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Chiral solitons in nuclei: Saturation, EMC effect and Drell-Yan experiments 
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Abstract 

The Chiral Quark-Soliton model of the nucleon contains a mechanism for an attractive interaction between nucleons. This, 
along with the exchange of vector mesons between nucleons, is used to compute the saturation properties of infinite nuclear 
matter. This provides a new way to asses the effects of the nuclear medium on a nucleon that includes valence and sea quarks. 
We show that the model simultaneously describes the nuclear EMC effect and the related Drell-Yan experiments. 



One frontier of strong interaction physics lies in the in- 
termediate range of length scales available to present ex- 
periments where neither the fundamental theory. Quan- 
tum Chromodynamics (QCD), nor its low energy effec- 
tive theory, Chiral Perturbation Theory, have useful per- 
turbative expansions. Neither fundamental quarks nor 
point-like hadrons provide a complete description, so in- 
cluding the non-perturbative information that hadrons 
arc bound states of valence quarks in a polarized vac- 
uum is necessary. One way to probe these intermedi- 
ate length scales and this non-perturbative physics is to 
examine the short distance structure of a large object. 
The prime example is the European Muon Collaboration 
(EMC) effect [l| where the short distance (~ 5 GeV, or 
^ 10~^ fm) structure of nuclei differs from that of a col- 
lection of free nucleons. This measurement showed that 
bound nucleons are different than free ones, and implied 
that the medium modifications could be significant for 
any nuclear observable 0] • Indeed, a recent paper 0] ob- 
tains evidence for a medium modification of the elastic 
proton form factor. 

Our central concern is the depletion of the nuclear 
structure function F2{x) in the valence quark regime 
0.3 % X 0.8. While the general interpretation is that 
a valence quark in a bound nucleon has less momen- 
tum than in a free one, corresponding to some increased 
length scale, the specific mechanism for this has eluded a 
universally accepted explanation for 20 years 0, 0, 0, ■ 
A popular explanation is the so-called 'binding' effect 
which originates from a possible mechanism in which 
mesons binding the nucleus carry momentum. An im- 
portant consequence is that the mesonic presence would 
enhance the anti-quark content of the nucleus IdilSl] • Such 
an effect has not been seen in Drell-Yan experiments 
in which a quark in a proton beam annihilates with an 
antiquark in a nuclear target producing a muon pair. 
Furthermore, relativistic treatments, including the light- 
cone approach needed to obtain the nucleon structure 
function, of the binding effect with structureless hadrons 
fail [ifll, luJ, , suggesting that modifications of the 

internal quark structure of the nucleon are required to 
explain the deep inelastic scattering data. 



Any description of the EMC effect must be consistent 
with the constraints set by both deep inelastic scatter- 
ing and Drell-Yan data. Thus a successful model must 
include antiquarks as well as quarks, and show how the 
medium modifies both the valence and sea quark dis- 
tributions. Our purpose is to provide a mechanism for 
that modification within the Chiral Quark-Soliton (CQS) 
model [l^ lisl Hil [l7| . This phenomenological model has 
many desirable qualities: the ability to describe a wide 
class of hadron observables with surprising accuracy, the 
inclusion of antiquarks, positivity of Generalized Parton 
Distributions, and a basis in QCD The model also 
predicted the recently discovered 0"^ exotic baryon 
resonance ^^HOl- Here we show how the model describes 
nuclear saturation properties, reproduces the EMC effect, 
and satisfies the bounds on nuclear antiquark enhance- 
ment provided by Drell-Yan experiments. 

The CQS model Lagrangian with (anti)quark fields 
'i/'jV'j ^i^d profile function 0(r) is 



(1) 



where 6(r oo) = and 0(0) = — tt to produce a soli- 
ton with unit winding number. The quark spectrum con- 
sists of a single bound state and a filled negative energy 
Dirac continuum; the vacuum is the filled negative con- 
tinuum with = 0. The wave functions in this spectrum 
provide the input for the quark and antiquark distribu- 
tions used to calculate the nucleon structure function. 

We work to leading order in the number of colors 
[Nc = 3), with Nf = 2, and in the chiral limit. While 
the former characterizes the primary source of theoretical 
error, one could systematically expand in Nc to calculate 
corrections. We take the constituent quark mass to be 
M = 420 MeV, which reproduces, for example, the A^-A 
mass splitting at higher order in the Nc expansion, and 
other observables |l7|. 

The theory contains divergences that must be regu- 
lated. We use a single Pauli-Villars subtraction as in 
Ref. liJl because we follow that work to calculate the 
quark distribution functions. The Pauli-Villars mass is 
determined by reproducing the measured value of the 
pion decay constant, /tt = 93 MeV, with the relevant di- 
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vergent loop integral regularized using Mpy — 580 MeV. 
This regularization also preserves the completeness of the 
quark states [2H . 

The nucleon mass is given by a sum of the energy of a 
single valence level {E^), and the regulated energy of the 
soliton {Eq^ equal to the energy in the negative Dirac 
continuum with the energy in the vacuum subtracted) 



Mjv = NcE"" + Eq{M) 



EeiM 



PV) 



The field equation for the profile function is 



Q{r) — arctan 



Pl{r) ' 



(2) 



(3) 



where pi and are the quark scalar and pseudoscalar 
densities, respectively. 

The dependence of nucleon properties on the nuclear 
medium is incorporated in the model by simply letting 
the quark scalar density in the field equation Q con- 
tain a (constant) contribution arising from other nuclc- 
ons present in symmetric nuclear matter. This models 
a scalar interaction via the exchange of multiple pairs of 
pions between nuclcons. We take the scalar density to 
consist of three terms: 1) the constant condensate value 
{'4'tp)o (in the vacuum or at large distances from a free 
nucleon), 2) the valence contribution and 3) the con- 
tribution from the medium which takes the form of the 
convolution of the nucleon and valence quark scalar 
densities as in the Quark-Meson Coupling (QMC) model 

[H Hi 113 

p1{r) ^ + p:{r) + / d\'p^{r ~ r')pUr'). (4) 



We take the pseudoscalar density to have only the valence 
term p^^ ~ p^^; the two other contributions analogous to 
the first and third terms of Eq. Q vanish due to sym- 
metries of the QCD vacuum and nuclear matter. These 
approximations to the densities neglect the precise form 
of the negative continuum wave functions in Eq. |(2Jl. The 
resulting free nucleon profile function has no discernible 
difference from a fully self-consistent treatment, demon- 
strating the excellence of this approximation. 

We take the vacuum value of the chiral condensate in 
Eq. (0J to be free parameter, but in the single Pauli- 
Villars regularization, the scalar and pseudoscalar densi- 
ties contain a divergence that cancels in the ratio Eq. Q 
[2^ . A finite value of the ratio is obtained by normal- 
izing the densities so that the cancellation occurs (by 
dividing the other terms in the numerator and denomi- 
nator by the same divergent quantity). This yields a free 
nucleon mass that is independent of {iptp)o (as is neces- 
sary because it is divergent in single Pauli-Villars regu- 
larization), and a medium contribution that enters only 
through the ratio p^ / {il!ip)o. While the vacuum value of 
the condensate does not vary by definition, the effective 
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FIG. 1: Binding energy per nucleon as a function of Fermi 
momentum fit to B = —15.75 MeV at the minimum for 
(iii>)o = -(225 MeV)^ (long dashed), -(210 MeV)=' (short 
dashed) and —(200 MeV)"^ (solid). The box and shaded region 
are the experimental uncertainty |2l| in the binding energy, 
density and compressibility of nuclear matter. 



condensate {■ipip)o + p^ S{kp), where S{kp) is integral of 
(see Eq. Q), falls ~ 30% at nuclear density. This 
is consistent with the model independent result |29l | that 
predicts a value 25-50% below the vacuum value. 

The nucleon scalar density is determined by solving 
the nuclear self-consistency equation 



= 4 



Mn{p, 



N\ 



(27r)3 ^P-f Mw(pf)2 



(5) 



The dependence of the nucleon mass, and any other prop- 
erties calculable in the model, on the Fermi momentum 
hp enters through Eq. (j^l. Thus there are two coupled 
self-consistency equations: one for the profile, Eq. lO, 
and one for the density, Eq. ij^l. These are iterated un- 
til the change in the nucleon mass Eq. (|2Jl is as small as 
desired for each value of the Fermi momentum. We use 
the Kahana-Ripka (KR) basis [2^, with momentum cut- 
off and box size extrapolated to infinity, to evaluate the 
energy eigenvalues and wave functions used as input for 
the densities, nucleon mass, and quark distributions. 

A phenomenological vector meson (mass m„ = 
770 MeV) exchanged between nucleons (but not quarks 
in the same nucleon), is introduced 12'^ to obtain the 
necessary short distance repulsion which stabilizes the 
nucleus. The resulting energy per nucleon is 



E 
A 



psikf 



d^k 
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(27r: 



3 ^k^ + MN{kpY + -^PB{kp). 



(6) 

We now present the results. The mass of a free nucleon 
is computed to be MN{kp^Q)^ 1209 MeV. The - 30% 
difference is as expected in the model at leading order in 
Nc- We evaluate the nucleon mass Eq. Q and energy 
per nucleon Eq. © as a function of hp for three values 
of the condensate. We plot B = E/A - Mn{Q) in Fig.Q 
where we choose the vector meson coupling to fit S = 

— 15.75 MeV at the minimum. 

1/3 

We use the value — ('0V')o = 200 MeV, and vector 
coupling (7^/47r = 10.55, which gives a Fermi momentum 
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^ in nuclear matter consistent with the 
= 1.35 ± 0.05 fm-i The compress- 



of kp = 1.38 fm 
known value kp 
ibility is if = 348.5 MeV which is above the experimen- 
tal value K = 210 ±30 MeV, but weU below the Walecka 
model value of 560 MeV. 

The isoscalar unpolarized distribution q{x) ~ u{x) + 
d{x) is the leading order term in Nc, with the isovec- 
tor unpolarized quark distribution u{x) — d{x) smaller 
by a factor ^ 1/A^c and set to zero. The distribu- 
tions are calculated using the KR basis at fc^ = and 
kp = 1.38 fm~^ almost exactly as in Ref. [i^ where the 
quark distribution is given by the matrix element 



n 

(7) 

with the regulated sum taken over occupied states. The 
eigenvalues are determined from diagonalizing the 
Hamiltonian, derived from the Lagrangian in the 
KR basis. The vector meson exchange is not explicit in 
Eq. (UI) because the initial, intermediate and final states 
of the struck and spectator quarks experience the same 
vector potential, as demanded by consistency. Thus we 
include the interaction of the debris of the struck nucleon 
with the residual nucleus The antiquark distribu- 

tion is given by q{x) = —q{—x) where the sum is over 
unoccupied states. We use the exact sea wave functions, 
and not the approximation used in Eq. The use of 
a finite basis causes the distributions to be discontinu- 
ous. These distributions are smooth functions of x in the 
limit of infinite momentum cutoff and box size, but nu- 
merical calculations are made at finite values and leave 
some residual roughness. This is overcome in Ref. 0| by 
introducing a smoothing function. We deviate from their 
procedure, and do not smooth the results; instead we find 
the subsequent one-loop perturbative QCD evolution [s^l 
to be sufficient. 

These distributions are used as input at a scale of Q = 
Mpv ^ 580 MeV for evolution to Q = 5 GeV in the 
case of the quark singlet distribution q^{x) = q{x) + 
q{x) oc F2 {x)/x at leading order in Nc- The EMC ratio 
function is defined to be 



R{x,Q^) 



F^{x,Q\kp) 
AF^{x,Q\kp = Qy 



(8) 



F^{x,Q\kp) = ^ dyf{y)F^{xly,Q\kp). 

The nucleon momentum distribution, following from a 
light-cone approach, for any mean field theory of nuclear 
matter llCl is 



fiv) 



{l + Ap-y)e{y-l + ^p) [A|, - (1 - yf] 



_ _ (9) 

where Af = kp/MN and Mn = Mjv(O) - 15.75 MeV. 
The antiquark distribution q{x) is evolved to Q = 




Aq 



FIG. 2: The EMC (top) and the Drell-Yan (bottom) ratios 
at scales Q = 5 GeV and 10 GeV, respectively, for nuclear 
matter. The data are for Iron (em pty boxes) and Gold (filled 
boxes) from SLAC-E139 (top) |3l|| for Iron (empty boxes) and 
Tungsten (filled boxes) from FNAL-E772 (bottom) 0. 



10 GeV for use in the Drell-Yan ratio q^/Aq, analogous 
to Eq. (P. The EMC and Drell-Yan ratios are plotted 
in Fig. 13 While the data shown in Fig. |2] are for large, 
but finite, nuclei, our calculation reproduces the trend of 
both sets of data. It falls shghtly below the SLAC-E139 
data due to the higher density of nuclear matter. 

In Fig. 131 we show the quark, antiquark, singlet and 
valence {q" = q — q) quark distributions weighted by x 
for a free and bound nucleon at a scale Q = 5 GeV. 
There is a large depletion in the bound nucleon valence 
distribution in Fig.|3 that, if used to calculate the EMC 
ratio (|H1), produces too large an effect. This large effect 
is comparable to that of the QMC model impulse ap- 
proximation calculation or the Guichon model 24] which 
only include valence quarks. This valence effect is miti- 
gated by a small enhancement in xq, consistent with the 
Drell-Yan data, so that the singlet distribution has only 
a moderate depletion consistent with the EMC effect. 

A simple picture in terms of the uncertainty principle 
is available. The influence of the nuclear medium on the 
nucleon causes the root mean square radius of the baryon 
density to increase by 2.4%. This corresponds to a de- 
creased momentum, and a depletion of the bound struc- 
ture function relative to the free one. This swelling is 
consistent with a < 6% increase as constrained by quasi- 
elastic inclusive electron- nucleus scattering data |32| | , and 
the recent polarization transfer measurement Q . 

We ignore the effects of shadowing, which occur when 
the virtual photon striking the nucleus fluctuates into a 
quark- antiquark pair over a distance ~ 1/2M]\[X exceed- 
ing the inter-nucleon separation. This causes a depletion 
in the structure function for x ^ 0.1 and is relatively well 
understood H, ^ |M 1^ and so we do not reiterate those 
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to calculate possible medium modifications of a variety 
of observables that could be measured in experiments. 

We would like to thank the USDOE for partial support 
of this work, and S. D. Ellis for useful comments. 
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FIG. 3: Clockwise from the top left the distributions 
xq{x),xq(x),xq^ {x),xq^ (x) in a free (dashed) and bound 
(solid) nucleon at a scale Q — 5 GeV. The valence distri- 
bution q^ {x)/2 is also shown in gray in the lower right graph. 



results. Additionally, we ignore contributions from quan- 
tum pion structure functions, which in this model propa- 
gate through constituent quark loops, and would modify 
the behavior at small x. These loops are suppressed by 
0{1/Nc), and are not treated at leading order. 

The present model provides a intuitive, qualitative 
treatment that maintains consistency with all of the free 
nucleon properties calculated by others 0, 0| . It gives 
a reasonable description of nuclear saturation properties, 
reproduces the EMC effect, and satisfies the constraints 
on the nuclear sea obtained from Drell-Yan experiments 
with only two free parameters: {'ipip)o and g^. 

The central mechanism to explain the EMC effect is 
that the nuclear medium provides an attractive scalar in- 
teraction that modifies the nucleon wave function. This 
is also the dominant mechanism in the QMC model ap- 
proach to the EMC effect 01 and also similar to the 
quark delocalization approach js^. The improvements 
given here are the explicit computation of the effects of 
the medium on the antiquark distributions so that consis- 
tency with the Drell-Yan data could be verified, and the 
reduction of the number of input parameters and model 
assumptions. Our extension of the Chiral Quark-Soliton 
model to nuclear matter provides a new, consistent way 
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